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Abstract 

We study models of quantum statistical mechanics which can be solved 
by the algebraic Bethe ansatz. The general method of calculation of 
correlation functions is based on the method of determinant repre- 
sentations. The auxiliary Fock space and auxiliary Bose fields are 
introduced in order to remove the two body scattering and represent 
correlation functions as a mean value of a determinant of a Fredholm 
integral operator; the representation has a simple form for large space 
and time separations. In this paper we explain how to calculate the 
mean value in the auxiliary Fock space of asymptotic expression of 
the Fredholm determinant. It is necessary for the evaluation of the 
asymptotic form of the physical correlation functions. 
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1 Introduction 



In this paper we use an example of the quantum nonlinear Schrodinger equa- 
tion, i.e. the one-dimensional Bose gas with delta-function interactions, in 
order to illustrate the development in the theory of quantum correlation 
functions. The correlation function of local fields in this model was studied 
in previous papers and its determinant representation was obtained 

in the paper QXJ| . The representation of the correlation function in terms 
of the Fredholm determinant of a linear integral operator is the basis of 
our approach. The differential equations for the Fredholm determinant had 
been obtained in M. They are directly related to the classical nonlinear 
Schrodinger equation. These differential equation were solved in the asymp- 
totic regime of large space and time separations; the simplified asymptotic 
form of Fredholm determinant was obtained in 0—0. This expression is 
an operator in an auxiliary Fock space. In order to find the asymptotics of 
the correlation function, one should calculate the vacuum mean value of this 
expression. It is a necessary step in the calculation of asymptotics of physical 
correlation functions. The problem of evaluation of the vacuum expectation 
(mean value) is a combinatorial problem, closely related to the procedure of 
the normal ordering in the quantum field theory. In the present paper we 
are studying just this problem. 

We remind in brief the basic definitions of the model under consideration 
for the reader's convenience. The quantum nonlinear Schrodinger equation 
can be described in terms of the canonical Bose fields ip(x,t) and i()'(x,t) 
(x € R) obeying the equal time commutation relations 



Here < c < oo is the coupling constant and h is the chemical potential. 
The spectrum of the model was first described by E. H. Lieb and W. Liniger 



i/}(x,t),i[)i(y,t) =8{x-y). 




The Hamiltonian and momentum of the model are 




and 
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The Lax representation for the corresponding classical equation of 

motion 

i—i> = [V, H] = -— </> + 2ctfiH> - hip, (1.4) 

was found by V. E. Zakharov and A. B. Shabat |7|]. The quantum in- 
verse scattering method for the model was formulated by L. D. Faddeev 
and E. K. Sklyanin §. 

The quantum nonlinear Schrodinger equation is equivalent to the Bose gas 
with delta-function interactions. In the sector with A" particles the Hamilto- 
nian of the Bose gas is given by 

N pfi 

^ = -E^ + 2c E S( Zj -z k )-Nh. (1.5) 

j=l aZ j l<j<k<N 

In this paper we shall consider the thermodynamics of this model. The 
partition function and the free energy of the model are defined by 

Z = tr e~~ = e~T. (1.6) 

The free energy F has been explicitly represented in terms of the Yang- Yang 
equation || 

^ = " 2 - h -tJ-j» sux-^ 1 + e ~^y (i7) 

F=~— ln(l + e t). (1.8) 
The correlation function, studied in this paper is defined by 

MO, 0W(x, t)) T = L JJ . (1.9) 

trie t) 

The plan of the paper is the following. In Section 2 we shall remind 
the reader of the asymptotic expression for the Fredholm determinant that 
represents the correlation function. We shall also describe its dependence on 
the quantum fields. Section 3 is devoted to the main results of the paper. In 
it we develop a technique for the evaluation of the mean values in auxiliary 
Fock space. It is related to the problems of the normal ordering in quantum 
field theory. In Section 4 we use this technique in order to calculate the mean 
value of the asymptotic expressions for the Fredholm determinant. 
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2 Asymptotics of the Fredholm Determinant 

In order to find the determinant representation of the correlation functions 
one has to introduce an auxiliary Fock space and three Bose fields ip(\), 
0(A), and $(A), which are linear combinations of annihilation and creation 
operators p(X) and g(A): 

^(A) = qty(A) +fty(A), 

0(A) = g^(A)+^(A), 
$(A) = g $ (A)+ P$ (A). 



The operators p(X) annihilate the Fock vacuum 

p(A)|0) = 0, (2.2) 
and the corresponding creation operators q(X) annihilate the dual vacuum 

(0|g(A) = o. (2.3) 

We shall also use the function 

h(X,fi) = (A- n + ic)/ic, (2.4) 
which enters into the non vanishing commutators 

MA), ?*(/*)] = -b^(A), qM\ = In (^j) . (2-5) 

[p%b(X),qM] = [p*(A),^(/x)] = tp^(A),g^(^)] = ln(/i(A,/i)/i(/i, A) J. (2.6) 



The relation of these quantum fields to the ones 0a 2 > 4>D\ used in the papers 
1-i are 

0(A) = A2 (A) - Dl (A), $(A) = A2 (A) + Dl (A). (2.7) 
The vacuum vector is normalized by unity (0 1 0) = 1. 



4 



The quantum fields ( [2.1Q are linear combinations of the three canonical 
Bose fields. The derivative of the field if)(\) also will be important: 







V'(A) = — ^(A) = ^(A)+ P ;(A). 



(2.8) 



Nonzero commutation relations between the derivatives of annihilation op 
erators p(A) and creation operators q(X) are: 

2ic 

[p^(A),g0(/x)] = MA),^,0)] 



A — /i + ic 



c 2 + (A-/i) 2 ' 

2 / 1 \ 2 

I + 



H — X + ic 



[p*(A),g^)] = [p^(A),^(/x)] = -[^(A),g$(/i)] = ^zjA^T^ ' 



(2.9) 
(2.10) 
(2.11) 



It is worth mentioning that quantum fields (|2.1| ) belong to the same 
Abelian sub-algebra. They all commute: 

= [^(A),0(/i)] = hKA), $(//)] = o, 

(2.12) 

[0(A), $(//)] = [0(A), 0(//)] = MX), $(//)] = 0. 

This property plays a very important role for the calculation of vacuum mean 
values in auxiliary Fock space. 

In the paper |T| the correlation function of local fields of the quantum 
nonlinear Schrodinger equation was represented as a mean value of a deter- 
minant of an integral operator, depending on the fields (|2.1| ). At large space 
x and time t separation the determinant simplifies to [ffl: 



(ii>{0,0)^(x,t)) T = (0\Q(x,t) 



1+0 u 



10), 



(2.13) 



where Q(x,t) is an operator in auxiliary Fock space 



Q(x,t) = C([0(A)],[$(A)])(2t) 



(u-l) 2 /2 ^(Xo)-itXl-iht 



1 f°° 

x exp <! — / dX (\x - 2Xt\ - i sign(A - X )ijj (X)) 

Z1l J-oo 



x In 



1 -tf(A) 1 + e 



3 0(A) sign(A-Ao) 



(2.14) 
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Here C is a smooth bounded functional, which may depend on x and t only 
through the ratio x/2t = Aq, which remains fixed. The Fermi weight $(A) is 
defined by 

(2.15) 



m = (i 



3 e(A)/T 



and 



v = J- ln{ [l - 0(A O )(1 + e^ (Ao) )l [l - 0(A O )(1 + e^ (Ao) )l } . (2.16) 
In this paper our main purpose is to evaluate the mean value of the right 



hand side of ( |2.14| ). Due to the relations ( j2.5| ) and ( |2.6| ), the creation and 
annihilation parts of the fields <fi and $ commute with each other. Thus 
the only nonzero contribution to the vacuum mean value is provided by 
normal ordering of expressions containing the field ip{X). It is easy to find 
the contribution of the factor e^ -* (see ( 2.14| )). Indeed let us move this 
exponent to the left 



MM) 



(0|e 



After this one can move e p ^ x °^ to the right using obvious relations: 

M A o),0(A o )] =0, 

b^(A ), H = o, 

= ( MO + In 



and 



e ^(Ao) $(//) = ($(^) +i n [/ i (A ^)% ) Ao)])e 

e P^(Ao)| ) = |Q). 



Pv( A o) 



(2.17) 

(2.18) 
(2.19) 

(2.20) 
(2.21) 
(2.22) 

(2.23) 



Thus we arrive at 

(0\Q(x,t)\0) = (0\C([<f>(\)}, [$(A)])(2t)^ 1 ) 2 / 2 e- l4A o-^ 

f 1 f°° / 2fA A 

x exp ^ — / dX (\x- 2Xt\ - i sign(A - A ) [^'(A) + " ' " 



2tt J- 



x In 



1 - 0(A) 1 + exp sign(A - A ) U{\) + In 



(A-A ) 2 + c 
h(X, X 



h(X ,X) 



|0). (2.24) 
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Here the functional C([0], [$]) can be obtained from the functional C([0], [$]) 
by shifting of the arguments of the last one according to rules ( j2.20| ) and 

Now the right hand side of ( |2.24|) can be written in the following form 



(0\Q(x,t)\0) = (0|exp 



dA^'(A)/(A|0(A)) F([0],[$])|O). (2.25) 



Here 



/(a|0(A)) 



sign (A - A c 
2m 



x In 



1 - i?( A) 1 + exp sign(A - A ) 0(A) + In 



h(X,Xo) 
h(X ,X) 



and 



F([0],[$])=C([0(A)],[$(A)])(2t) 



(v-l) 2 /2 e -it\l-iht 



x exp < — / dX (\x — 2Xt\ — 2i sign(A — A c 
[ 2n J-oo \ 



X-Xr 



x In 



1 - (A) ( 1 + exp [sign(A - A ) (0(A) + In 



(A - A ) 2 + c 
h(X,X 



h(X , A) 

In the next section we shall evaluate the right hand side of (|2.25| 



(2.26) 



. (2.27) 



3 Evaluation of the Mean Value 

The main purpose of this section is the evaluation of the mean value 

(0|exp( r rfA^(A)/(A|0(A))|F([0],[$])|O). (3.1) 



Here complex function / becomes an operator, because it depends on quan- 
tum field 0(A). It is worth mentioned that the particular case of Q3.1|) , when 



/(A|0(A)) is a linear function of the field 0(A), was first considered in [ID|. 
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We remind the reader of the definitions 

^(A)=^(A)+ g ;(A), 

0(A)=p (A) + g (A), (3.2) 
$(A) = P$ (A) + g$(A). 

As usual, the relations p(X) |0) = and (0|g(A) = for all p and q are satisfied 
as well as the commutation relations 

b4( A )'^(X>] = £0^) = M A )>^0")]> ( 3 - 3 ) 

[p^(A), g$(/i)] = f (A, /i) = [p$(A), gj(/x)], (3.4) 
[^(A),^0)] =r)(\,/i). (3.5) 
The complex functions £(A,/i) and £(A,/x) are equal to 

e(A^) = ^A) = c2 + ( 2 ; C _^ )2 , ^^)- (A 2 _ A /x ^ c2 i (3-6) 

and 2 

»l(A, = T x * ) +f x 1 ■ ) • (3-7) 
\X — H + ic ) \X — fi — ic) 

However we do not use explicit expressions ( |3.6| ) and ( |3.7| ) in this section. 

We evaluate the mean value (|3.1|) in three steps. First, we consider an 
auxiliary problem — the scalar case. 



3.1 Scalar Case 

Consider the following mean value 

(o|e {^(A)/(^))} F([0U$])|o) _ m 

Here F([(j)], [$]) is a smooth functional, depending on the fields and $. A 
complex function f(z) becomes an operator- valued function while its argu- 
ment 0(/i) is an operator. As usual such an expression should be understood 
as a formal Taylor series; therefore, we assume that f(z) is holomorphic 
within some circle \z\ < po- 
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The parameter a is a complex number such that the following restriction 
holds: 

II P° 
\z\ = p < 



a 



/(z£(A, n) 



< z 



for 



ie(A,A*)i 



(3.9) 



The complex numbers A and fi are fixed. 

Let us decompose the exponent in a Taylor series, 



3 {<#'(A)/(^))} = y a™/ 



E^(^(A))>(W). 



(3.10) 



n=0 



We would like to emphasize that here we essentially used the commutativity 
of the fields ip' and (ft. 

In order to calculate (0| (if}'(\)) n we use the Cauchy integral representation 

p zV'(A) 

(3.11) 



(V>'(A)) r 



dz' 



nl 



z=0 



2m 



dz- 



yU+l 



\ Z \=P 



The evaluation of (0|e 2 ^'^ A ' ) is a standard problem in quantum field theory 

(Oje^'W = e z2r? ( A ' A )/ 2 (0|e^ (A) . (3.12) 



After substituting this into (|3.11|) , we obtain 

p z 2 r,{\,\)/2 



(0i(V'(A))" 



2vri 



-(0|e 



(3.13) 



Further substitution into ( |3.8|) and ( |3.10|) gives the expression for the mean 
value 



00 a n 



E 

?1=0 



2vri 



dz 



e z 2 V (X,\)/2 



(O|e^( A V n (0(/x))F([0],[$])|O). (3.14) 



y n+l 



2|=p 



Recall that creation and annihilation parts of the fields <fi and $ commute 
with each other; therefore, we have 

$\^+ {X) nm)m\, mm = Mt&wrbMmfo], mio) 



r(^(A,/*))(O|e^WF([0],[$])|O). (3.15) 
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If we substitute this into ( p.!4|) then we may calculate the sum with respect 
to n. Due to the restriction ( |3.9|) this series is absolutely convergent: 



(O|e^'W /( ^»F([0],[$])|O) = — 

\z\=p 



,2 

dz 



e z 2 V (X,\)/2 



z - atf(z£{\,iiyj 



x(O|e^ (A) F([0],[$])|O). (3.16) 
Due to the relation ( |3.9| ) and the Rouche theorem, the equation 

z- a/ =0 (3.17) 

has exactly one zero of the first order z = z = z (a) in the domain \z\ < 



p. Let us emphasize that (|3.17| ) is a classical equation where only complex 
functions are involved (i.e., no quantum operators). Thus, taking the integral 
with respect to z, we arrive at 



(0|e<«~>» F([A m)|0) = ^(0|^.'^([y])|0) (3 18) 

The right-hand side may be further simplified 

e^ {A ' A)/2 (O|e z < (A) F([0], [$])|0) = (O|e^' (A) F([0], [$])|0). (3.19) 



Let us recall once more that in the original formula ( |3.8| ) /(</>(//)) was a 
quantum operator since it depends on the quantum field 0(/x). The result of 
the above calculation shows that this function may be replaced by a complex 
number z . 

It is worth mentioning, that one can analytically continue the result ob- 
tained with respect to a into the domain where the inequality (|3.9| ) is not 
valid. We propose the theorem, 

Theorem 3.1 

(0{e {^M} mi mo) = W^WMM (3 . 20) 

1 - af'[z ^(X,id)j 
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Here the complex number zq can be found from the equation 

z = af(zo£(\,[i)), 

and 



(3.21) 



(3.22) 



Remark. The equation ( [3.21 ) may have many solutions if we do not 
impose the restriction ( |3.9| ). In this case one should choose the solution 
z = z (a), with the property: 



z (a J 



0. 



a=0 



(3.23) 



Let us remind the reader that no operators are involved in equation ( ft.21| ) 
since f(z) is a complex function. This is complex equation for the complex 
number zq. 



3.2 Matrix Case 



The method of calculation of the mean value described above can be easily 
generalized for more complicated case. Namely let us consider the example: 



JV 



(0| exp £ V'(A fc )AMA fc )) F([0], [<&]) |0). 



(3.24) 



-fc=i 



It is clear that we may find the mean value in an analogous fashion to that 
for a scalar case. Let us briefly describe the main steps of the corresponding 
derivation. 

First, we have 



N 



N oo 



exp \J2^'(X k )fMX k )) = II E -Tf/'i A,))"'/;' SoiA,)). (3.25) 



T) ■ ' 
j=ln J = 3 ' 



The normal ordering of ip'i^j) niay be performed in a similar fashion to ( |3.11| ) 
and (gT2l) : 



N 

in 



OO 1 



n,=o n r 



(2ni 



N dz- 

rtj+1 
j=l Zj 
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x exp J WniXj, A fc ) 1 (0|e^^i^ (Afc) . (3.26) 

[ 1 i,k=i J 

Here each of the integrals is taken by circle, which appears to be a common 
domain of analyticity of the functions fk(z). 
Next we find the mean value 

N N N N 

( | e £^ J] F([4>], [$])|0) = n /r(E^(A„ A fc ) 



fe=l 



fc=l i=l 



X 



( O | e Sf=i^( A '=)F([0],[$])|O). 



(3.27) 



Now we substitute this into the expression for our mean value (p.24|) and 
sum with respect to each n^. 



N 



X 



e 1 /2E^=i^ z ^' A '=)(o|e^^i^ (Afc) F([0],[$])|O) 



nJLi 



/j ( Em=l ^m£(Am, Aj) 



(3.28) 



In order to take the integral we introduce the complex numbers z® as 
solutions of the system 



N 



m=l 



We also define the matrix M: 
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Mj k = 5jk — -z — fk( £ m £(A m , Aa 
m =i 



(3.29) 



(3.30) 



After evaluation of the z« integration, we obtain the result 



(0| exp £ ^(A*)/*(0(A*)) F([0], [$])|0) 



-fc=i 



exp{Eti^(A fc )^°}F([0],[$])|O) 



detM 



(3.31) 
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As in scalar case, the system ( |3.29| ) may have many solutions. In this 
case one can consider the replacement f k — > af k . After this, it is necessary 
to choose the solution of ( ^.29[ ), which approaches zero as a — > and to 
continue this solution to the point a — 1. 

3.3 Continuous Case 

In order to evaluate the mean value (|3.1| ) let us consider the continuous limit 
of flOIQ : 

\ k+1 = A fc + A, (3.32) 

f k (<f>(\ k )) = Af(X k \(f>(X k )). (3.33) 

4 = Az(A fc ), (3.34) 

and take the limit A — > 0. We notice that the system of constraints in ( |3.29| ) 
turns into the integral equation: 

Z (X) = f(\\ J^dfi z{ti)£(ji,\)). (3.35) 

Furthermore, the matrix Mj k becomes an integral operator with the kernel 

M(X,ii) = S(X- l i)-j^rf(X\ jTtte z(s)£(s,X)). (3.36) 
Equation (|3.31| ) has the following continuous limit: 
(0|exp{£^A ^(A)/(A|(0(ai)))f(M,[*])|O) 

= (detM)- 1 (0|exp<{ / dXip'(X)z(X) \ F([<f>), [$])|0). (3.37) 



Equation (|3.35|) -( |3.37|) are the main result of this section. We would like to 
emphasize that the commutativity of the quantum fields ( |2.1|) is extremely 
important for obtaining of this result. 

The following calculations are trivial. We have 

(0| exp ( r dX^'{X)z{X) 



eXP {^/°° dXdfX ^( A '^^ A ) Z ^)} (°l ex p{/°° dAp^(A)z(A)J . (3.38) 
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The action of the operator exp{p'^} on the functional F leads to the shift of 
the arguments of the last one (see ( j2.20|) , ( |2.21|) ). We find 



(0| exp dXp'^X)z(X)} F([<j>], [$])|0) 



d\z{X)i(X,n) 



(3.39) 



In the next section we shall use these results in order to evaluate the mean 
value of the asymptotic expression ( 2.24J) . 



4 The Mean Value of the Leading Term 

In order to find the mean value of the leading term of asymptotics (|2.24j) , 
we need only to substitute the concrete expressions ( |2.26| ), fl2.27Q into the 
equations (|3^35| ) - (|3~37|) and ( ^39l) . 

An integral equation for the z-function is 



z{X) = -— sign(A - A ) In {1 - (A)X(A, A )} 

Z7T 



(4.1) 



where 



X(X, A ) = 1 + exp jsign(A - A ) ( In ^ Ao) + 



h(X ,X) 



2icz{n) 



c 2 + (A - fi) 2 



The integral operator M has the kernel 



(4.2) 



M(A, /i) = 8(X -/i)- ^sign(A - A )-^- In {1 - #(X)X(X, A )} . (4.3) 

After evaluation of the variation derivative one should set z(X) equal to the 
solution of (O). 



The functional -F([0], [$]) is given in ( p. 27] ). It consists of three factors: an 
exponential factor, a power law correction with respect to t, and a constant 
factor, which depends only on the ratio x/2t = A . Using ( |3.39| ) and the 
integral equation (|4.1| ), we find the exponential factor to be 



exp 



dX {x - 2Xt)z{X) 



(4.4) 
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(4.5) 



The next factor is a power law correction 

where the expression for v follows from ( |2.16|) : 

v = ±- In { [l - #(A )(1 + e-^ (Ao) )l [l - #(A )(1 + e*< A °>)l } . (4.6) 

Instead of </>(A ) we substitute the integral expression, 

2ic 



0(A O ) -> tt(A ) 



r*(A), 



(4.7) 



C 2 + (A -A) 2 ' 

where the function z(X) is a solution of ( |4.1| ). So the power law correction 
becomes 

(2t)(-D 2 /2 ) 

where 



i? = J- In { [l - i?(A )(l + e^ A °))l [l - #(A )(1 + e u(Ao) )l } . 



2tt 

Finally, the constant term is equal to 

f 1 /"OO 

^ = C([i*0*)],[t;0i)])e^e 3 q) 
where 



(4.8) 
(4.9) 



rfA 



2 J — OO 

2(A-/i) 



cZAcZyu, r/(A, fjL)z(X)z(fi) 
z(X). 



-oo c 2 + (A - A) 2 
Thus we obtain for the mean value (|2.24j ) of the operator Q(x,t) 



(4.10) 
(4.11) 



(0\Q{x,t)\Q)=g-e 



(det M)~ l 



x exp < — z 



dA (z - 2At)*(A) 



(4.12) 



This is the main result of the paper. 

Expression ( |4.12j ) is the leading term of the asymptotic evaluation of the 
Fredholm determinant representing the correlation function Q. In the next 
publication we shall study the corrections. The mean value of corrections can 
contribute to the leading term of the asymptotic behavior of the correlation 
function. 
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